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O Abstract 

$_! The elastic radial deformation of multiwall carbon nanotubes (MWNTs) under hydrostatic pressure is investigated within the con- 
CIh tinuum elastic approximation. The thin-shell theory, with accurate elastic constants and interwall couplings, allows us to estimate 
the critical pressure above which the original circular cross-section transforms into radially corrugated ones. The emphasis is 
placed on the rigorous formulation of the van der Waals interaction between adjacent walls, which we analyze using two different 
approaches. Possible consequences of the radial corrugation in the physical properties of pressurized MWNTs are also discussed. 
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1. Introduction 

Nanomaterials are tiny platforms on which a beautiful in- 
terplay between structure and property can be appreciated. 
For the last decade, advanced nanofabrication techniques have 
put out various nanostructured materials with novel geometry 
HI S B EL many of which exhibit unprecedented proper- 
ties not seen in macroscopic structures. Among such nanos- 
tructures, carbon nanotubes have drawn great deal of attention. 
The salient feature of carbon nanotubes is their mechanical ro- 
bustness and resilience |@]: Due to the extremely large stiff- 
ness, for instance, their thermal conductivity becomes higher 
than even that of diamond n. Besides, their flexibility in 
bending twisting ||9j, [10|, LLil, [12|], radial compression 

|[T3ill4ll5Lll6l[l7ll . and the associated variations in the physical 
properties hold promise for developing nanoelectromechanical 
devices 1 18, T^]. Towards successful implementations of such 
ideas, computational studies have been playing a vital role in 
complementing experimental observations, often difficult and 
incomplete for nanomaterials |20| . 

In view of structure-property relations, radial deformation is 
expected to elicit untouched behavior of multiwall carbon nan- 
otubes (MWNTs); see Fig.[T]for a microscopic image of a typ- 
ically synthesized MWNT 12111 . The radial deformation dis- 
turbs both the equal spacings between the concentric walls in 
MWNTs and the in-plane hexagonal lattice within each wall. 
In particular, spatial modulation in the wall-wall separation may 
enhance (or hinder) the electron charge transfer between neigh- 
boring walls if they locally get close to (or pull away from) 
each other 11221 12311 : this implies significant changes in various 
observables of MWNTs under external pressure. Despite the 
broad interest, theoretical efforts on the mechanics of radially 
pressurized MWNTs remain limited 124l I25L 12611 . This scarcity 
is mainly because atomistic-based simulations of MWNTs hav- 
ing many walls are computationally very expensive. 



In the present work, we employ a powerful alternative, the 
continuum elastic thin-shell theory approach, to analyze the sta- 
ble cross-sectional shapes of MWNTs under hydrostatic pres- 
sure. Such approaches are often plagued with inaccurate or in- 
consistent elastic moduli. Here, we carefully select the intra- 
wall and wall-wall elastic constants for an accurate prediction 
of the critical pressure p c above which the original circular 
cross-section transforms into radially corrugated ones. In the 
corrugation modes, each wall exhibits a wavy structure in the 
circumferential direction along the tube axis. Such the pressure- 
induced corrugation is attributed to the mismatch in the me- 
chanical stability between the flexible outmost walls with large 
tube diameters and rigid innermost walls with small diameters. 



2. Methods 

2.1. Mechanical energy 

Although atomistic simulations may provide precise estima- 
tions of physical quantities in general, they often demand huge 
computational resources in systems of interest. On this back- 
ground, the thin shell theory based analysis for the carbon 



nanotube mechanics has long been developed II 1 3L 1271 l28l 
30l 31]. Along the continuum thin-shell method, a MWNT is 
mapped onto a set of N continuum elastic hollow tubes of radii 
r, : (1 < i < N). A point on the circle corresponding to the cross- 
section of the zth tube is described by (x, y) = (r,- cos 8, r,- sin 8) 
in terms of the polar coordinates; 8 is the circumferential angle 
around the tube axis. Underpressure p, the point moves to 

(x*,y*) = ( [n + Ui(e,p)] cos e-Vi(6,p) sin0, (1) 
[r t + Ui(6, p)] sin 8 + Vi(G, p) cos 8 ) . (2) 

If the deformation amplitudes, and Vj, are sufficiently 
small, the mechanical energy of the ith tube per unit axial length 
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Figure 1 : Transmission electron microscopy image of a MWNT obtained via 
chemical vapor decomposition. The number of concentric walls, N, and the 
innermost tube diameter, D, are estimated as 12 and 3.0 nm, respectively. After 
Ref. Q. 



is given by 



U 



(0 



Here, e, and /c, are respectively the in-plane and bending strains 
of the ;th wall, both depending on u„ v, and their derivatives 
with respect to 9 13211 . The constant C denotes the in-plane 
stiffness, D the fiexural rigidity, and v the Poisson ratio of each 
wall. 

For quantitative discussions, the values of C and D must 
be carefully determined. In the conventional thin-shell theory 
for macroscopic objects, C and D are related to the Young's 
modulus E of the wall and its thickness h as C - Eh, D = 
£7i 3 /[12(l - v 2 )]. However, for carbon nanotubes, the wall is 
made out of a monoatomic graphitic layer and consequently the 
notion of a wall thickness becomes elusive. Hence, the macro- 
scopic relations for C and D noted above fail since there is no 
unique way of defining the thickness of the graphene wall 133 1. 
Thus the values of C and D should be evaluated ab-initio from 
direct measurements or computations of carbon sheets, with- 
out reference to the macroscopic relations. In actual calcula- 
tions, we substitute C = 345 nN/nm, D = 0.238 nN- nm, and 
v = 0.149 along with the prior work 028I1 based on the density 
functional theory. It should be noted that the values of C and 
D are essentially tube-diameter dependent. Neverthelss, such 
the dependences become negligible when the tube diameter ex- 
ceeds 1 nm, above which the elastic constants of carbon nan- 
otubes converge to those of a planar graphene sheet |28]. On 
this background, we will take into account only the nanotubes 
whose diameters are larger than 1 nm, which allows to fix the 
values of C and D as noted above. 

The stable cross-sections of MWNTs under p minimize the 
mechanical energy U of the whole system that is described by 
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with u' = du/dO. The term IT™ accounts for the van der Waals 
(vdW) interaction energy of adjacent pairs of walls, and O is 
the negative of the work done by p during cross-sectional de- 
formation. Note that U is a functional of m, and v,; therefore, 
the variational method allows to obtain the optimal solutions of 
iij and v, that minimize U under a given p. 

For N » 1, the outside walls have large diameters, and con- 
sequently are very flexible and susceptible to mechanical insta- 
bilities under radial pressure. The contrast in the radial rigidity 
between the outermost walls and the innermost ones triggers 
a non-trivial cross-sectional deformation observed in radially 
pressurized MWNTs. In fact, the linearized buckling analysis 
leads to the conclusion that immediately above a critical pres- 
sure p c , the circular cross section of MWNTs becomes radially 
deformed as described by 134 ] 



Uj{0) = uf\p c ) + 5pi(n) cos n8, v,(f?) = 6vj(n) sin nd. (7) 

The solution O represents a wavy structure of a MWNT's 
cross-section, called the radial corrugation with a mode index 
77. The integer n indicates the wave number of the corrugated 
walls, being uniq uely determined by the one-to-one relation be- 
tween n and p c 13411 . We will find below that n depends sys- 
tematically on and the innermost tube diameter D = r\ as 
presented in Fig. [3] 

2.2. Wall-wall interaction coefficient 

It is known that in the axial buckling of MWNTs, the buck- 
ling load is sensitive to the vdW interaction between adjacent 
walls l35ll . Likewise, p c is thought to depend on the strength of 
vdW interactions; this is the reason why we address the rigorous 
formulation of the interaction. Using both a discrete-pairwise 
summation and a continuum approach, we derive the coeffi- 
cients c,j in Eq. (0 through a first order Taylor approximation 
of the vdW pressure i3~5l[3al associated with the vdW potential 
V(d) = 4e[(cr/rf) 12 - (cr/c/) 6 ]. Here, d is the distance between a 
pair of carbon atoms, 2^ 6 cr = 0.3833 nm is the equilibrium dis- 
tance between two interacting atoms, and s = 2.39 meV is the 
well depth of the potential 0711 . The resulting equilibrium spac- 
ing between neighboring walls, used here to define the geome- 
try of the MWNCTs, is 0.3415 nm. The derivative F = dV/dd 
is the force between the two atoms. 

There exist several continuum models for the vdW interac- 
tions lf35ll36il3~7ll . In Ref. ll35Tl . expressions for the pressure and 
the vdW interaction coefficients were obtained by integrating 
the continuum vdW force and its derivative on curved wall sur- 
face, while disregarding the vectorial nature of the force. The 
significance of the vectorial nature of the force was addressed 
in Ref. ll36ll . where analytical expressions for the vdW pressure 
were obtained by considering only the component of the vdW 
force normal to the wall. It was emphasized in Ref. (36] that 
for a two-walled carbon nanotube, the pressure exerted on inner 
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Figure 2: Tube-radius dependence of the derivatives of pressure (top) and the 
vdW interaction coefficients c,-j (bottom). 



wall is different from the pressure on the outer wall. The pres- 
sures on the inner and outer walls for a concentric two-walled 
tube with radii r mn and r out are given below (with positive signs 
for compression): 



'"out r ! Inn r 

Pmn = a — /_ and p out = a /+, 

r inn '"out 



(8) 



where a = 3nscrp c 2 /32 with p c = 38.18 nm 2 being the area 
density of carbon atoms. We used the notations: 

f ± = 231/3 11 {jE n +E U )- 160/? 5 (yE 7 ± E 5 ) , (9) 

where f3 = cr/(r out + r inn ), y = h/(r om + r inn ), h = r out - r inn , 
E >» = IT^ 1 - £ 2 sin 2 c?r ffl/2 <f6) and jfc = 4r mn r out /(r mn + r out ) 2 . 

In the following, we obtain analytical expressions for cy lin- 
earizing the formula ([8]) for the pressure. The formula takes 
into account correctly the normal-to-wall component of vdW 
forces, and avoids common assumptions like p; nn = p oat Ir29ll 
and pinn^inn = Pout^out B 3 811 ; therefore, the resulting c,j will be 
also free from unnecessary assumptions. 

In the present work, infinitesimal deformation is considered. 
Hence, the linearzed pressure is needed for vdW energy calcu- 
lation. Note that the vdW energy depends quadratically on the 
change in spacing between two adjacent walls. Consider two 
consecutive walls with radii r, and r, + i, where the subscripts 
i and ; + 1 correspond to inn and out respectively. The vdW 
energy stored due to symmetric perturbation along the positive 
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Figure 3: (color online) N-dependence of critical pressure curve p c . Imme- 
diately above p c , the original circular cross section of MWNTs gets radially 
corrugated. 



direction of pressure is given by 

r m r ln ( Ah Ah\ 

J„ (-^t-^t)* (10) 

where r m = (r, + r, + i)/2 is the mean radius and is the 
vdW pressures on the ith wall. The corresponding linearized 
pressure is given by dpij+i/dh\ . In Eq. dTOb . r m d0 describes 
the length of the infinitesimal element on which the pressure 
is acting. Using the linearized pressure and comparing with 
Eq. ©, we get the expressions for vdW coefficients as 



cu+i — 7 



dp i dp i+ 



dh 



dh 
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See Appendix for the derivatives of df ± /dh. Note that the c ;j is 
symmetric. 

These analytical expressions for the stiffness coefficients as- 
sociated with the vdW interactions obtained with the continuum 
approximation are compared with the coefficients obtained us- 
ing a discrete pairwise calculation. In the later approach, the 
pressure coefficients are computed by summing the component 
normal to the wall of the linearized vdW forces acting on an 
atom a, i.e. the normal component of 2/>(<9f 7 dd)(dd / dh)\^ a b ), 
where the linearized force dF/dd is evaluated at the actual dis- 
tance between the pair of interacting atoms a and b, d a -\,. In 
Ref. |j35]] > the linearization is done at the equilibrium distance of 
the potential. Noting that most of the interacting atom pairs are 
not at the equilibrium distance, it is apparent that the lineariza- 
tion should be done at the actual distance between each pair of 
interacting atoms. The discrete pressure can be obtained by di- 
viding the force exerted on one atom by the area per atom. The 
vdW force on one atom is obtained by considering the inter- 
actions between this atom with all the atoms on a neighboring 
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Figure 4: (color online) Independence of normalized critical pressure curve p c . 
Inset: Cubic behavior of the p c curve when N = 1 . The dashed line showing 
the power-law p c oc D~ 3 is a visual reference. 



wall. The effect of the relative lattice placement of interacting 
walls is averaged out by integrating the interactions over a sym- 
metric triangle of the representative hexagonal cell of graphene. 

The normalized pressure derivatives (dpjj/dh)/(&7Tso-p c 2 ) 
and normalized vdW interaction coefficients Cjj/(8nso-p c 2 ) are 
plotted against normalized radius {r mn + r out )/cr in Fig. [2] A 
good match is observed between the discrete-pairwise and the 
continuum approach. The interaction coefficient for two hori- 
zontally placed graphene walls, cf ap em , j s obtained from the 
second derivative of the analytic continuum interaction poten- 
tial [37]. For larger radii, c,j approaches to c ^ pheae 7 corrobo- 
rating our results for nanotubes. 



3. Results and discussion 

Figure [3]plots p c as a function of N for various values of D. 
An initial increase in p c at small N (except for D — 1 .0 nm) 
is attributed to the enhancement of radial stiffness of the en- 
tire MWNT by encapsulation. This stiffening effect disappears 
with further increase in N, resulting in decay or convergence 
of p c {N). It is noteworthy that MWNTs practically synthesized 
often show D larger than those presented in Fig. [3] In fact, the 
MWNT depicted in the image of Fig.[T]gives D = 3.0 nm, for 
which p c (N) lies at several hundreds of MPa as estimated from 
Fig. [3] Such degree of pressure applied to MWNTs is easily ac- 
cessible in high-pressure experiments, supporting the feasibility 
of our theoretical predictions. The D-dependence of p c for vari- 
ous fixed N is provided in Fig.|4]in a normalized manner; When 
N - 1, we obtain p c oc D 1-3 in a gree ment with previous studies 
on single- walled nanotubes [Q~3,l39]. It also deserves comment 
that a radial pressure large enough to cause corrugation can be 
achieved by electron-beam irradiation; the self-healing nature 
of eroded carbon walls gives rise to a spontaneous contraction 
that exerts a high pressure on the inner walls 

Figure [5] shows the stepwise increases in the corrugation 
mode index n, defined by Eq. (7). For all D, the deformation 
mode observed just above p c abruptly increases from n — 2 to 



Figure 5: (color online) Stepwise increase in the corrugation mode index n. 



n > 4 at a certain value of N, followed by the successive emer- 
gence of higher corrugation modes with larger n. These suc- 
cessive transitions in n at N » 1 originate from the mismatch 
in the radial stiffness of the innermost and outermost walls. A 
large discrepancy in the radial stiffness of the inner and outer 
walls results in a maldistribution of the deformation amplitudes 
of concentric walls interacting via vdW forces, which conse- 
quently produces an abrupt change in the observed deformation 
mode at a certain value of N. 

A possible physical consequence of radial corrugation is a 
pressure-driven change in quantum transport of carriers in ra- 
dially corrugated nanotubes. It is known that mobile carriers 
whose motion is confined to a thin curved layer behave differ- 
ently from those on a conventional flat plane because of the 
curvature-induced electromagnetic field [43, 44, 14511 . Another 
interesting issue is the effect of atomic lattice registry on non- 
linear deformation {i.e. large deformation amplitude) regimes 
11461 l47l 14811 . In the latter case, the degree of commensurance 
in atomic structures between neighboring carbon layers plays 
a prominent role in determining the optimal morphology of 
highly corrugated MWNTs, in which the interwall spacings par- 
tially vanish. Crumpling or twisting caused by radial pressure is 
expected, and in fact were observed in preliminary calculations 
B 



4. Conclusion 



The thin-shell theory has been employed to derive the critical 
hydrostatic pressure p c above which the cross section of pres- 
surized MWNTs transforms into radially corrugated patterns. 
The value of p c lies on the order of several hundreds of MPa for 
practically synthesized MWNTs with Af » 1 and D > 3 nm, 
indicating the experimental feasibility of the prediction. The 
enegetically favored corrugation pattern strongly depends on N 
and D, where a larger N and smaller D give a larger corruga- 
tion mode index n. Our results provide a clue in developing 
MWNT-based devices operating under high pressure. 



4 



Acknowledgement 

We acknowledge M. Rahimi, E. Osawa, Y. Suda, K. Yakubo 
and T. Mikami for helpful discussions. HS is thankful for the 
financial supports provided by KAKENHI and Hokkaido Gas 
Co., Ltd. SG acknowledges the support of the Spanish Ministry 
of Science and Innovation through the Juan de la Cierva pro- 
gram. MA acknowledges the support of the European Research 
Council (FP7/2007-2013)/ERC grant agreement nr 240487 and 
the prize "ICREA Academia" funded by the Generalitat de 
Catalunya. Numerical simulations were partly carried out using 
the facilities of the Supercomputer Center, University of Tokyo. 

Appendix A. Functional forms of df±/dh 
It follows from Eq. (|9]l that 



dh 



= 23i/*»(2*i3+y^±^) 



(A.l) 



All E m functions obey the recurrence relation (m - 2)(1 - 
k 2 )E m = (m - 3)(2 - K 2 )E,„-2 - (m - 4)E m _4 from complete 
elliptic integrals of the first and second kind, K(k) and E(k), re- 
spectively; Ei = K(k) and £_i = E(k) in the current notation. 
Hence, the required derivatives of E m in Eq. ( IA.U are given 
below: 



dEs 
Ok 

0E-, _ 
ak \5(k 



dk 



dEr 



, 1 -3 (-3E - 7k 2 E + 2k 4 E + 3K- 2k 2 K - k 4 K] , 
1 -r ( 15 E + 58 k 2 E - 33 k 4 E + 8 k b E 

l) k V 

- 15 K + 2 k 2 K + + 17 k 4 K - 4 k 6 K) , 
1 n6 (-315 kE - 2309 k 2 E + 261 1 k A E 

315(-l+/t 2 ) k V 

-1S99 k 6 E + 760 k*E -128 k w E 

+3 15 kK + 419 k 2 K - 1403 k A K + 993 k 6 K 

-3SS k & K + 64 k w K), 

tt- (-3465 E - 31907 k 2 E + 44991 k 4 E 

Yk V 



ak 3465(-l+/t 2 )'/t 

-43633k 6 E + 26206 k*E - 8832 £ 10 £ + 1280 k 12 E 
+3465 K + 7652 k 2 K - 24562 k 4 K + 23204 k 6 K 
-13615 k*K + 4496 k w K - 640 k n K) . 

(A.2) 

Other derivatives are given by dk/dh\,- m = -hj(r 2 Ji) and 
dyldh\ Tm = l/r m . 
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